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Abstract Dendritic spines are small protrusions on a neuronal dendrite that are the 
main locus of excitatory synaptic connections. Although their geometry is variable 
over time and along the dendrite, they typically consist of a relatively large head 
connected to the dendritic shaft by a narrow cylindrical neck. The surface of the head 
is connected smoothly by a funnel or non-smoothly to the narrow neck, whose end 
absorbs the particles at the dendrite. We demonstrate here how the geometry of the 
neuronal spine can control diffusion and ultimately synaptic processes. We show that 
the mean residence time of a Brownian particle, such as an ion or molecule inside 
the spine, and of a receptor on its membrane, prior to absorption at the dendritic shaft 
depends strongly on the curvature of the connection of the spine head to the neck 
and on the neck's length. The analytical results solve the narrow escape problem for 
domains with long narrow necks. 



1 Introduction 



Recognized more than 100 hundreds years ago by Ramon y Cajal, dendritic spines 
are small terminal protrusions on neuronal dendrites, and are considered to be the 
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Fig. 1 Upper row (from left to 
right): The dendritic spine head 
is connected smoothly to the 
neck (the postsynaptic density is 
marked red) and (right) the 
connection is not smooth (from 
[6]). Lower row: 

A mathematical idealization of a 
cross section: A cross section of 
a sharp and non- sharp 
connection approximating the 
spine morphology. 





main locus of excitatory synaptic connections. The general spine geometry consists 
of a relatively narrow cylindrical neck connected to a bulky head (the round part in 
Figure 1). Their geometrical shape correlates with their physiological function [1-6]. 
More than three decades ago, the spine-dendrite communication associated with the 
particle transfer was already anticipated [7] to be mediated not only by pure diffu- 
sion but it was hypothesized to involve other mechanisms such as twitching. This 
idea was reinforced by the findings [8] that inside the spine, the cytoplasmic actin is 
organized in filaments, involved in various forms of experimentally induced synap- 
tic plasticity by changing the shape or volume of the pre- and postsynaptic side and 
by retracting and sprouting synapses. The fast dendritic spine contraction was finally 
confirmed in cultured hippocampal neurons [9] and consequences were studied theo- 
retically in [10-12]. Interestingly, a serial electron microscopy and three-dimensional 
reconstructions of dendritic spines from Purkinje spiny branchlets of normal adult 
rats allowed to relate spine geometry to synaptic efficacy [1]. This image reconstruc- 
tion approach leads to the conclusion that the cerebellar spine necks are unlikely to 
reduce transfer of synaptic charge by more than 5-20%, even if their smooth endo- 
plasmic reticulum were to completely block passage of current through the portion of 
the neck that it occupies. The constricted spine neck diameter was proposed to isolate 
metabolic events in the vicinity of activated synapses by reducing diffusion to neigh- 
boring synapses, without significantly influencing the transfer of synaptic charge to 
the postsynaptic dendrite [1]. 

Change of spine morphology can be induced by synaptic potentiation protocols 
[13-15] and indeed intracellular signaling such as calcium released from stores alters 
the morphology of dendritic spines in cultured hippocampal neurons. These changes 
in geometry can affect the spine-dendrite communication. One of the first quantita- 
tive assessment of geometry was obtained by a direct measurement [16] of diffusion 
though the spine neck. Using photobleaching and photorelease of fluorescein dextran, 
by generating concentration gradients between spines and shafts in rat CA1 pyrami- 
dal neurons, the time course of re-equilibration was well approximated by a single 
exponential decay, with a time constant in the range of 20 to 100 ms. The role of the 
spine neck was further investigated using flash photolysis of caged calcium [3, 17] 
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and theory [18], and the main conclusion was that geometrical changes in the spine 
neck such as the length or the radius are key modulator of the spine-dendrite commu- 
nication [12, 19, 20], affecting calcium dynamics. However, in all these studies, the 
nature of the connection between the neck and head was not considered. The theoreti- 
cal studies [19, 21] considered non- smooth connection only of the head to the narrow 
cylindrical neck (Figure 1) and did not account for any effect of curvature. This is 
precisely the goal of this article to investigate the consequences of this connection. 

The connection between the head and the neck is not only relevant for the three- 
dimensional diffusion, but also essential to the analysis of other synaptic proper- 
ties. Indeed, synaptic transmission and plasticity involve the trafficking of receptors 
[22-27] such as AMPA or NMDA receptors (AMPARs or NMDARs) that mediate 
the glutamatergic-induced synaptic current. Single particle approaches have further 
[28, 29] revealed the heterogeneity of two-dimensional trajectories occurring on the 
neuron surface, suggesting that there are several biophysical processes involved in 
regulating the receptor motion. In addition, the number and type of receptors that 
shape the synaptic current [23] could be regulated by the spine geometry. This ques- 
tion was further explored theoretically [30, 31], using asymptotic expressions for the 
residence time and experimentally [32] by monitoring the movements of AMPARs 
on the surface of mature neurons using FRAP. Employing a combination of confo- 
cal microscopy and photobleaching techniques in living hippocampal CA1 pyramidal 
neurons, a correlation between spine shape parameters and the diffusion and compart- 
mentalization of membrane-associated proteins was recently confirmed [33]. Lateral 
diffusion seems to be a constitutive process of AMPAR trafficking; it depends on 
spine morphology and is restricted by the spine geometry [34] . 

In this article, we develop a method for computing the NET from composite spine- 
like structures that consist of a relatively large compartment (head) £2\ and a narrow 
cylindrical neck of cross section \dQ a and length L (see Figure 2). Our connection 
formula is given as 



The connection between the two parts in the context of the NET problem was at- 
tempted in [21, 35] for the oversimplified model of a discontinuous connection. Here, 
we study a large class of connections and reveal the role of curvature in the spine- 
neck connection in regulating diffusion flux through narrow necks. More specifically, 
we study here the residence time of a Brownian particle from the spine head to the ab- 
sorbing end of the spine neck moving either on the surface or inside the spine. We use 
the results of [36, 37] for the mean first passage time (MFPT) to an absorbing bound- 
ary at the end of a cusp-shaped protrusion in the head. They account for the effects 
of curvature generated by the neck-head connection in the spine. The reciprocal of 
the MFPT is the rate of arrival (probability flux) of Brownian particles from the head 
to the dendrite [38]. We calculate the narrow escape time (NET) from spine-shaped 
domains with heads connected smoothly and non- smoothly to the neck. 




L 1 \Qi\L 
2D + \3Q a \D 



(1) 
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Fig. 2 The composite domain 
consists of a head Q\ connected 
by an interface 3 £2; to narrow 
neck ^2 • The entire boundary is 
reflecting, except for a small 
absorbing part dQ a . 




2 The NET from a domain with a bottleneck 

We consider two- and three-dimensional composite domains £"2 that consist of a head 
£2i connected through a small interface 3 £2; to a narrow cylindrical neck £22. The 
boundary of £2 is assumed reflecting to Brownian particles, except the far end of £22, 
denoted 3£2 fl , which is absorbing. For example, in Figure 2, the interface 3 £2/ is the 
black segment AB and the absorbing boundary dQ a is the segment CD at the bot- 
tom of the strip. The NET from such a composite domain cannot be calculated by 
the methods of [39-42], because the contribution of the singular part of Neumann's 
function to the MFPT in a composite domain with a funnel or another bottleneck is 
not necessarily dominant. The method of matched asymptotic expansions requires 
different boundary or internal layers at a cusp-like absorbing window than at an ab- 
sorbing window which is cut from a smooth reflecting boundary (see [43-46]). The 
methods used in [21, 35] for constructing the MFPT in a composite domain of the 
type shown in Figure Id are made precise here and the new method extends to the 
domains of the type shown in Figure lc. 

First, we recount some basic facts about the NET [35, 39-41, 43-45, 47, 48]. The 
NET is the MFPT of a Brownian trajectory to a small absorbing part of the bound- 
ary of a domain, whose remaining boundary reflects Brownian trajectories. Refined 
asymptotic formulas for the NET were derived in [42, 46, 49, 50], and were used to 
estimated chemical reaction rates. 

Consider Brownian motion x(t) in a sufficiently regular bounded domain £2, 
whose boundary 3 £2 consists of a reflecting part 3 £2 r and an absorbing part 3 £2^ . 
The expected lifetime of x(t) in £2, given x(0) = x e £2, is the MFPT v(x) of x(t) 
from x to 3£2 fl is the solution of the mixed boundary value problem [38] 

Av(x) = forxe£2 (2) 

D 

v(x) = 0 forxedQa (3) 

dv(x) 

— — =0 forxe3£2 r , (4) 
dn 

where dv(x)/dn is the normal derivative at the boundary point x. If the size of the 
absorbing part 3 £2^ of the boundary is much smaller than the reflecting part 3£2 r , the 
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MFPT r = v(x) is to leading order independent of x e Q a and can be represented by 
the Neumann function N(x, y) as 

1 C f dv(x) 

r = v(y) = -- / N(x,y)dx- / N(x,y)— — dS x . (5) 

The sum of the integrals is independent of y e £2 a outside a boundary layer near 3 Q a . 
The Neumann function is a solution of the boundary value problem 



A x N(x, y) = -8(x - y) + — for x, y e Q 



SN(x,y) 
dn x 



= 0 forxedQ,yeQ 



(6) 
(7) 



and is defined up to an additive constant [39, 47]. 
2. 1 The MFPT from the head to the interface 

In the two-dimensional case considered in [40] the interface 9 £2; is an absorbing 
window cut from the smooth reflecting boundary of Q\, as in Figure Id. The MFPT 
?x^dQi is the NET from the reflecting domain Q\ to the small interface 8Qi (of 
length a), such that s = 7r|9£2i|/|9£2i | = 7ta/\dQ\ | <^ 1 (this corrects the definition 
in [40]). It is given by 



IQiL la^il 

ln- 



•0(1) 



7tD 7t\dQi\ 

for x e £2\ outside a boundary layer near 8Qi . 
In particular, if Q\ is a disk of radius R, then for x = the center of the disk, 



(8) 



R^ 
D 



R 1 

log - +21og2+ - + 0(e) 
a 4 



averaging with respect to a uniform distribution of x the disk 



?x^dQi 



R 
~D 



2 r 



R 



1 



log- +21og2+- + 0(e) 



(9) 



(10) 



When the interface dQi (of length a) is located at an algebraic cusp with radius of 
curvature R c (see Figures lc and 2), the MFPT is given in [36, 37] as 



r = 



4Dj2a/R c 



(1 + 0(1)) fore«|3^|. 



(11) 



In the case of Brownian motion on a spherical head of the surface of revolution ob- 
tained by rotating the curve in Figure Id about its axis, is a sphere of radius R 
centered at the origin, connected to Q2 by a circle dQi centered on the north-south 
axis near the south pole, with small radius a = R sin 8/2. The domain ^2 is a right 
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Fig. 3 A surface of revolution 
with a funnel. The z-axis points 
down. 
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cylinder of radius a connected io Q\ at dQi, and the absorbing boundary dQ a is the 
circle of radius a at the bottom of the cylinder. The MFPT from Q\ to dQi is given 
in [37,41,51-53] as 



2R 2 
~D 



Tx^dQt =— log— -f , (12) 



where 0 is the angle between x and the south-north axis of the sphere. 

A surface of revolution generated by rotating a curve about an axis, as in Figure 3, 
with a funnel of diameter s can be represented parametrically as 

x=r(z)cosO, y = r(z)sin(0), 0 < 0 < 2tt, 0 < z < A, (13) 

where the axis of symmetry is the z-axis with z = 0 at the top of the surface and 
z = A at the end of the funnel, r is distance to the z-axis, and r =r(z) is the equation 
of the generating curve. We have 

r{z) = 0{^Jz) nearz = 0 

r( z ) =a + r a (A- z) l+a (\ + o(\)) for a > 0 near z = A, (14) 

where I has dimension of length. For a = 1 the parameter i is related to the radius of 
curvature R c at z = A by i = 2R C . For a > 0 [37] 

2£> (l+^sin^ 

where \Q\ \ is the entire area of the surface. In particular, for a = 1 we get the MFPT 

l^il 

Tx^dQi ~ ——7===. (16) 

4D^2a/R c 

The case a = 0 corresponds to a circular cap of a small radius a cut from a closed 
surface. 

The MFPT of Brownian motion from a solid ball Q\ to a disk 3^- of small radius 
a near the south pole is given by [42] 



l^il 
4aD 

^ Springer 



(2 7? / a R 
1 + — log - + o - log - 



(17) 
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(note that the MFPT is independent of x up to second order, see [46]). For a general 
three-dimensional domain, Q\ the MFPT to a circular cap dQi cut from a smooth 
boundary is given by [42] 



IQll 

4a D 



1 + 



L x + Rx 
lit 



TV 



1/2 



log 



\8Qi 



la^il 



log, 



la^il 



+ 



0(1) 



(18) 



where L x , are the principal curvatures at a point x, and | = na 2 is the area 
of the circular cap. 

When the interface 3 Qi is a circular disk of radius a at the end of an axisymmetric 
solid funnel, the MFPT is drastically affected and changes to 



1 (R c 



3/2 



R C D 



(1+0(1)) fora«i? c , 



(19) 



where R c is the radius of curvature at the end of the funnel [37]. 



3 Connecting a head to a narrow neck 

We consider Brownian motion in a domain £"2 that consists of a head, which is a reg- 
ular bounded domain Q\, and a narrow neck Q2, which is a right circular or planar 
cylinder of length L, perpendicular to the boundary 8Q\, and of radius a (see Fig- 
ure 2). Thus, the interface dQi between the head and the neck is a line segment, a 
circle, or a circular disk, depending on the dimension. We assume that dQ\ — dQi 
is reflecting and that the other basis of the neck, dQ a c 3^2, is absorbing for the 
Brownian motion. The length (or area) \dQi \ is given by 



|3^| = 



a for a line segment 

2ica for a circle (20) 
it a 1 for a disk. 



The MFPT r x ^ d Q a can be represented as [54], [38, Lemma 10.3.1, p. 388] 

*X^dQ a = tx^dQi + tdQi^dQa, (21) 

where the MFPT TsQ^dQa is ?x^dn a , averaged over dQi with respect to the flux 
density of Brownian trajectories in Q\ into an absorbing boundary at dQi (see [38] 
for further details). 

First, we calculate rdQ.^sQ a an d the absorption flux at the interface. In the nar- 
row neck ^2 the boundary value problem (2)-(4) can be approximated by the one- 
dimensional boundary value problem 

Du zz = -1 
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w(0)=0, u(L) = uh, 
where the value at the interface u(L) = uh is yet unknown. The solution is given by 

z 2 

u(z) = -— + Bz, (22) 

so that 

L 2 

u(L) = u H = -— + BL, (23) 

which relates the unknown constants B and u h • The constant B is found by multi- 
plying Equation 2 by the Neumann function N(x,y), integrating over Q\, applying 
Green's formula, and using the boundary conditions (3) and (4). Specifically, we ob- 
tain for all y e dQi 

1 f f dv(x) 1 f 

u(j) = --/ N(x,y)dx- N(x,y)-^dS x + — v(x)dx. (24) 

Approximating, as we may, v(y) ~ u(L) and using (23), we obtain 

1 f f dv(x) 

BL = / N(x,y)dx- N(x,y)—-dS 

1 f 



2D 



(25) 



+ / v(x)dx . 



Because v(x) is the solution of the boundary value problem (2)-(4) in the entire 
domain Q = Q\ U ^2, the meaning of (25) is the connecting rule (21), where 

?x^dn a = 77^-7 f v(x)dx (26) 
l"i| Jq x 

tdn^dQa =u(L) (27) 

1 r r dv(x) 

f x ^ d Q i =-- / N(x,y)dx- / N(x 9 y)— —dS x . (28) 

Equation 26 gives the MFPT, averaged over Q\. The averaging is a valid approxi- 
mation, because the MFPT to dQi is constant to begin with (except in a negligible 
boundary layer). Equation 27 is the MFPT from the interface to the absorbing end 
dQ a of the strip, and (28) follows from (5). 

Matching the solutions in Q\ and Q2 continuously across dQi, we obtain the total 



flux on dQi as 



f 3u(jk 
= D — — 



dv(x) 

' dS x = -(\Qi\ + \Q 2 \), (29) 
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Noting that dv(x)/dn = -u f (0) = -B, we get from (20) and (29) that 

for a line segment 



B = - 



m l 

aD D 



Ifill 




L 


2naD 


D 


mil 




L 




D 



for a circle 

for a circular disk. 



(30) 



Finally, we put (21)-(30) together to obtain 



\&i\L 



(31) 



2D \8Q a \D 
The MFPT Tx-^bq. is given by (8)-(19) above. 

3.1 The NET from two- and three-dimensional domains with bottlenecks 

The expression (31) for the NET from a domain with a bottleneck in the form of 
a one-dimensional neck, such as a dendritic spine, can be summarized as follows. 
Consider a domain £2 with head £"2i and a narrow cylindrical neck ^2 of length L 
and radius a. The radius of curvature at the bottleneck in smooth connecting funnel 
is R c . In the two-dimensional case 



l^ii |ani| o(i) l 2 \Qi\L 

In I I I 

ixD a D 2D aD 

planar spine connected to the neck at a right angle 
TrlGi! []& . _ L 2 \Qi\L 



— ^(l+ 0 (l)) + + 
D V a 2D 2itaD 

planar spine with a smooth connecting funnel 



1^1 



I02 



sin : 



IfiilL 



(32) 



2ixD sinf 2D 2naD 
spherical spine connected to the neck at a right angle 

2D 2tt<2£> 



2Z) (1 + or) sin 



1 +a 



spherical spine with a smooth connecting funnel, 



where is the radius of the sphere, a = R sin 5/2, and 0 is the initial elevation angle 
on the sphere. If \Q\ | ^> aL and L ^> a, the last term in (32) is dominant, which is 
the manifestation of the many returns of Brownian motion from the neck to the head 
prior to absorption at dQ a (see an estimate in [19]). The last line of (32) agrees with 
the explicit calculation in [37]. 

The NET of a Brownian motion from a three-dimensional domain £2 with a bot- 
tleneck in the form of a narrow circular cylinder of cross- sectional area na 2 is given 
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Fig. 4 Left: The NET of Brownian motion on a sphere with a bottleneck connected by a smooth funnel to 
the neck (dashed line), and with a non-smooth connection (continuous line). Right: The NET of Brownian 
motion in a ball with a bottleneck connected by a smooth funnel to the neck (dashed line), and with a 
non-smooth connection (continuous line). 



by 



>d£2 a 



\Ql\_ 

4a D 



a i R 
1 + — log- 
ix R a 



0(1) L 1 \Qi\L 

J J L 1 _ 

D 2D na 2 D 



solid spherical head of radius R connected to the 
neck at a right angle 

1/2 




(33) 



D 2D na 2 D 
a general head connected to the neck at a right angle 



J_(*c 

V2\a 



— — (i + + — + — 

R C D 2D 7ta z D 



a general head connected smoothly to the neck by 
a funnel, 



where R c is the curvature at the cusp. The order 1 term can be computed for the 
sphere using the explicit expression of the Neumann-Green function [46]. Figures 4 
and 5 show the NET for various parameters, such as the neck length and radius. 

Finally, the influence of the neck length on the residence time is shown in Figure 5: 
changing the neck length modulates dramatically the residence time. Interestingly, the 
geometry of the connection affects much significantly the dimension two rather than 
that the three dimensional Brownian motion. 
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1 




0 I , ' ' 0' 

0.5 1 1.5 2 0.5 1 1.5 2 

Neck length L(fim) Neck length L(/im) 

Fig. 5 Effect of modulating the neck length. The geometry is the same as in Figure 4. 



4 Discussion and conclusion 

We have shown that the mean residence time (or flux) of Brownian particles inside a 
spine-like structure or on its surface depends strongly on the geometrical properties 
of both head and neck. Surprisingly, it also depends strongly on the smoothness of 
the connection between the two. 

The application to a freely diffusing AMPA receptor, which is responsible for the 
excitatory synaptic current, shows that its motion on the spine membrane is strongly 
restricted by dendritic spine geometry. Our results can be used to estimate the res- 
idence time of the receptor on the membrane if interactions with any scaffolding 
molecules are neglected (the latter are mostly concentrated in a local microdomain 
called the PSD). Using Equation 32 for non-smooth geometry (Figure lb) and for 
a spherical head of radius R = 1 /xm, a neck length L = 1 /xm, a neck radius 
a = 0.1 /xm, and a diffusion coefficient D = 0.1 /xm 2 /^, we obtain from the third 
line of Equation 32 that the residence time is T2d ~ 260 s, while for a smooth con- 
necting geometry (Figure la) with a curvature radius of R c = 1 /xm at the connecting 
neck-head T2d~1150s (line 2 in Equation 32). Evidently, the residence time is more 
than doubled, leading to the conclusion that there is a significant difference between 
the function of spines with smoothly and non- smoothly connected necks. We con- 
clude from this analysis that an AMPA receptors that do not interact with the PSD 
stay on a typical dendritic spine between one and a half to 5 min on average and this 
residence time is controlled mostly by the geometrical properties of the spine. 

We now consider the residence time for freely diffusing particles such as 
molecules, mRNA, and ions (e.g., calcium) inside a dendritic spine. For a calcium 
ion, the diffusion constant is about D ca = 400 /xm 2 /s [12]. Calcium ions that exit the 
spine only to the dendrite shaft at the end of the neck, but not through exchangers, 
give the following residence time estimates. Using formula 33 for a non-smooth con- 
nection between the spine head and the neck, we obtain that t^d ~ 195 ms (line 1), 
while for a smooth connection with radius of curvature 1 /xm, the residence time is 
T3D ~ 820 ms (line 3 in Equation 33). Interestingly, the mean residence time is tripled 
from the non- smooth to a smooth connection. 
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A remaining open question is to extend the present analysis to the case where many 
binding molecules can trap receptors. This effect should be expected to significantly 
increase the residence time inside a dendritic spine, as has already been observed in 
[55] for the case of a receptor inside the PSD. The present mathematical analysis of 
the residence time provides a solution to the narrow escape problem for domains with 
bottlenecks [21, 35]. Other generalizations of this study is to include the dynamics of 
many receptors [30, 56] or/and to study dendritic trafficking [57]. 

There are many other factors that affect the spine-dendrite communication with re- 
spect to calcium. This includes calcium pumps, endogenous buffers, calcium stores, 
the number and rates of exchangers. These mechanisms affect the residence of cal- 
cium in spines [58-60] and it would be interesting to add them in the present analysis. 
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